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ABSTRACT 

This manuscript is an updated version of Kalogera et al. (2004) published in ApJ Letters to correct 
our calculation of the Galactic DNS coalescence rate. The details of the original erratum submitted 
to ApJ Letters are given in page 6 of this manuscript. 

We report on the newly increased event rates due to the recent discovery of the highly relativistic 
binary pulsar J0737-3039 (Burgay et al. 2003). Using a rigorous statistical method, we present 
the calculations reported by Burgay et al. , which produce a coalescence rate for Galactic double 
neutron star (DNS) systems that is higher by a factor of 5 — 7 compared to estimates made prior to 
the new discovery. Our method takes into account known pulsar-survey selection effects and biases 
due to small-number statistics. This rate increase has dramatic implications for gravitational wave 
detectors. For the initial Laser Interferometer Gravitational- Wave Observatory (LIGO) detectors, the 
most probable detection rates for DNS in-spirals are one event per 10 — 630 yr; at 95% confidence, we 
obtain rates up to one per 3 yr. For the advanced LIGO detectors, the most probable rates are 10 — 500 
events per year. These predictions, for the first time, bring the expectations for DNS detections by the 
initial LIGO detectors to the astrophysically relevant regime. We also use our models to predict that 
the large-scale Parkes Multibeam pulsar survey with acceleration searches could detect an average of 
four binary pulsars similar to those known at present. 

Subject headings: binaries: close - stars: neutron - gravitational waves - methods: statistical 
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1. INTRODUCTION 

For almost 30 years, close double neutron star (DNS) 
binary systems have been known to exist in the Galaxy 
as a small subset of the observed radio pulsar popula- 
tion (Hulse & Taylor 1975; Wolszczan 1991). These sys- 
tems lose orbital energy due to the emission of gravita- 
tional waves (Taylor & Weisberg 1989, 2003; Stairs et 
al. 1998); the associated orbital in-spiral continues un- 
til the binary system coalesces, resulting in a burst of 
gravitational waves. DNS in-spirals are prime targets for 
gravitational-wave detection by the ground-based inter- 
ferometers Laser Interferometer Gravitational- Wave Ob- 
servatory (LIGO; Abramovici et al. 1992), GEO (Danz- 
mann et al. 1995), and VIRGO (Caron et al. 1997). 
Event rate estimates are very important for the devel- 
opment of gravitational-wave interferometers (Thorne & 
Cutler 2002). They are based on estimates of Galactic 
rates and their extrapolation throughout a survey vol- 
ume (Finn 2001), given the source strength and instru- 
ment sensitivity. For DNS binaries, Galactic rate esti- 
mates have been obtained using two very different meth- 
ods. One is purely theoretical and involves models of bi- 
nary evolution calibrated usually to the observationally 
determined supernova rate for the Galaxy. The other, 
more empirical, approach is based on the physical prop- 



erties of the close DNS binaries known in the Galactic 
field and modeling of radio pulsar survey selection ef- 
fects. For a review and details of both these approaches, 
see Kalogera et al. (2001, hereafter KNST) and refer- 
ences therein. The empirical method has generally pro- 
vided us with better constraints on the coalescence rate 
(KNST), although the uncertainty still exceeds two or- 
ders of magnitude. This is primarily due to (1) the very 
small number (only two until recently) of close DNSs 
known in the Galactic field with merger times shorter 
than a Hubble time and (2) the implicit assumption that 
this small sample is a good representation of the total 
Galactic population (KNST). 

Two recent developments make it appropriate to revisit 
the DNS merger rate calculations. First, the discovery 
of the 2.4 hr DNS binary PSR J0737-3039 in a large- 
area survey using the Parkes radio telescope (Burgay et 
al. 2003) brings the number of known DNS systems to 
merge in the Galactic field to three. With an orbital pe- 
riod of only 2.4 hr, J0737-3039 will coalesce in only 85 
Myr, a factor of 3.5 shorter than the merger time of PSR 
B1913+16. This immediately hints towards a possible 
significant increase of the coalescence rate (Burgay et al. 
2003). Second, a novel statistical method has been devel- 
oped by Kim, Kalogera & Lorimer (2003, hereafter KKL) 
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that automatically takes into account statistical biases 
inherent in small-number samples, like the relativistic 
DNS binaries, and in addition allows us to quantify our 
expectation that the actual DNS binary coalescence rate 
has a particular value, given the current observations. 

In this Letter, we use this statistical method and inves- 
tigate in detail the effect of this new DNS discovery on 
the estimates of Galactic DNS in-spiral rates and its im- 
plications for gravitational wave detection in this decade. 
We summarize the method in §[3 and discuss the result- 
ing Galactic in-spiral rate in §|3| In §^we use our models 
to make predictions for the expected number of DNS bi- 
naries that the Parkes Multibeam (PMB) survey (e.g. 
Manchester et al. 2001) could detect when acceleration 
searches are completed, and in §|Slwe discuss the impli- 
cations of our results for the detection rates of upcoming 
gravitational-wave detectors. 

2. METHOD FOR RATE CALCULATION 

Until recently, estimates of DNS coalescence rates pro- 
vided a range of possible values without any information 
on the likelihood of these values. KKL presented a newly 
developed statistical analysis that allows the calculation 
of a probability distribution for rate estimates and the 
determination of confidence intervals associated with the 
rate estimates. The method can be applied to any radio 
pulsar population (see KKL for the two close DNS known 
at the time and Kim et al. 2004 for close binaries with 
a pulsar and a massive white dwarf companion). Here 
we update the results of KKL, taking into account the 
recent discovery of the new DNS binary PSR J0737-3039 
(Burgay et al. 2003). 

The method is described in detail in KKL, but we 
briefly summarize the main elements here. The method 
involves the simulation of selection effects inherent in all 
relevant radio pulsar surveys and a Bayesian statistical 
analysis for the probability distribution of the in-spiral 
rate estimates. The small-number bias and the effect of 
the faint end of the pulsar luminosity function, previously 
identified as the main sources of uncertainty in rate esti- 
mates (KNST), are implicitly included in this analysis. 

For a model Galactic pulsar population with an as- 
sumed spatial and luminosity distribution, we determine 
the fraction of the total population which are actually 
detectable by current large-scale pulsar surveys. In or- 
der to do this, we calculate the effective signal-to-noise 
ratio for each model pulsar in each survey, and compare 
this with the corresponding detection threshold. Only 
those pulsars which are nominally above the threshold 
count as detectable. After performing this process on 
the entire model pulsar population of size A'tot, we are 
left with a sample of A"obs pulsars that are nominally de- 
tectable by the surveys. By repeating this process many 
times, we can determine the probability distribution of 
A^obs, which we then use to constrain the population and 
with a Bayesian analysis derive the probability expecta- 
tion that the actual Galactic DNS in-spiral rate takes on 
a particular value, given the observations. More details 
are given in § 2 of Lorimer et al. (1993) and in KKL. 

When this method was first developed (KKL), it was 
shown that, although the shape of the probability dis- 
tribution of rate estimates is very robust, the rate value 
at peak probability systematically depends primarily on 
the characteristics of the radio pulsar luminosity func- 



tion: its slope and the physical minimum luminosity of 
pulsars. Both of these are constrained by the general 
pulsar population (see Cordes & Chernoff 1997), but we 
explore the dependence of our results on the assumed 
values. 

Here we consider the same set of pulsar population 
models as in KKL, but we choose model 6 as our reference 
model in view of the recent discovery of very faint pulsars 
(for a review, see Camilo 2003). With the addition of 
the new DNS binary PSR J0737-3039, our calculations 
differ from those in KKL in two main ways: (1) the latest 
Parkes survey that led to the discovery of the new system 
(Burgay et al. 2003) is included, and (2) we calculate 
and account for the effects of Doppler smearing for DNS 
binaries akin to J0737-3039 by creating fake time series 
for a variety of orbital phases (see KKL for details) . Even 
for a ~ 4.5— min integration (Burgay et al. 2003), this 
effect alone reduces the average signal-to-noise ratio of a 
2.4-hr DNS binary by 35%. 

The statistical analysis presented in KKL has been 
extended to account for three systems (see Kim et al. 
2004). In our calculations we adopt a total lifetime for 
J0737-3039, defined as the sum of the current age and 
the remaining lifetime until the final coalescence, equal 
to 100 -f 85 = 185 Myr (Burgay et al. 2003). In the ab- 
sence of detailed beam observations for the new binary, 
we adopt a beaming factor of ~ 6 equal to the average of 
the two observationally constrained beams of the previ- 
ously known DNS binaries (KNST; see also Burgay et al. 
2003). We do note, however, that studies of other known 
recycled pulsars (the majority of them strongly recycled, 
spinning faster compared to the DNSs by about an order 
of magnitude) have shown that beaming fractions can 
vary significantly (Kramer et al. 1998). It is important 
to keep in mind that any uncertainties in the beaming 
factor proportionally affect the rate estimates, but not 
the rate increase factors derived here. 

3. GALACTIC IN-SPIRAL RATE 

For our reference pulsar model (with a radio luminos- 
ity function consistent with current pulsar observations; 
Cordes & Chernoff 1997; Camilo 2003), we find the most 
likely value of the total coalescence rate to be 7^ = 83 
Myr~^. The ranges of values at 68% and 95% confi- 
dence intervals are 40-170 and 20-290 Myr~^, respec- 
tively. The width of these ranges are somewhat smaller 
than previous estimates (the ratio between upper and 
lower limits at 68% confidence interval is 4.4 cf. 5.7 found 
by KKL), confirming the expectation that a bigger ob- 
served sample would reduce the uncertainty in the rate 
estimates (KNST). The new value for 7?, is a factor of 
6.4 higher than found by KKL. From the resulting prob- 
ability distribution shown in Figure ^ it is clear that 
J0737-3039 dominates the total rate over the other two 
systems. This is due to two separate factors: (1) The es- 
timated total number of DNS binaries similar to J0737- 
3039 (1700) is far higher than those of each of the other 
two systems (700 for B1913-f 16 and 500 for B1534-I-12). 
This is mainly due to the shorter pulsar spin and binary 
orbital period of J0737-3039, which results in a signif- 
icant Doppler smearing and efficiently "hides" them in 
the Galaxy (2) The total lifetime of J0737-3039 (185 
Myr) is significantly shorter than those of the other two 
(365 Myr for B1913+16 and 2.9 Gyr for B1534-)-12). 



IN-SPIRAL RATES OF DOUBLE NEUTRON STARS 



3 





'm 
C 

d) 
Q 

t -5 
'2 

(0 

o 



-10 



00 

o 



-15 



-20 



I I I I I III 1 1 I I I Mil 



Detection rate for initial LIGO (yr ') 
0.0001 0.001 0.01 0.1 1 

— I — I I I I I i[j 1 — I I I I iii| 1 — I I I 1 1 ii| 

B 1534+ 12 

- 



J0737-3039 




Table 1. Estimates for Galactic in-spiral rates and predicted LIGO 

DETECTION RATES (AT 95% CONFIDENCE) FOR DIFFERENT POPULATION MODELS 
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Galactic merger rate (Myr"') 

Fig. 1. — Probability density function that represents our ex- 
pectation that the actual DNS binary merger rate in the Galaxy 
(bottom axis) and the predicted initial LIGO rate {top axis) take 
on particular values, given the observations. The curves shown 
are calculated assuming our reference model parameters (see text). 
The solid line shows the total probability density along with those 
obtained for each of the three binary systems ( dashed lines) . Inset: 
Total probability density, and corresponding 68%, 95%, and 99% 
confidence limits, shown in a linear scale. 



We now explore our results for all other models con- 
sidered in KKL. Our main results are shown in Table 
1, where we have included a subset of models that re- 
flect the widest variations of the rates (as shown in KKL, 
variations in the space distribution of pulsars are not im- 
portant) . The main conclusions that can be easily drawn 
are: (1) The increase factor on the in-spiral rate is highly 
robust against all systematic variations of the assumed 
pulsar models and is strongly constrained in the range 
5-7; this is consistent with but somewhat lower than 
the simple estimate presented in Burgay et al. (2003). 
(2) The shape of the rate probability distribution also 
remains robust, but the rate value at peak probability 
depends on the model assumptions in the same way as 
described in detail in KKL (see Figs. 5-7 in KKL). 

4. PREDICTIONS FOR FUTURE DISCOVERIES 

As already mentioned, long integration times combined 
with very short binary orbital periods strongly select 
against the discovery of new binary pulsars. Specifically, 
in the large-scale PMB survey (e.g. Manchester et al. 
2001) with an integration time of 35min, the signal-to- 
noise ratio is severely reduced by Doppler smearing due 
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to the pulsars' orbital motion. Acceleration searches in 
the current re-analysis of the PMB survey (Faulkner et 
al. 2003) should significantly improve the detection effi- 
ciency to DNS binaries. 

Following Kalogera, Kim & Lorimer (2003), we calcu- 
late the probability distribution that represents our ex- 
pectation that the actual number of DNS pulsars with 
merger times shorter than a Hubble time (Nohs) that 
could be detected with the PMB survey takes on a partic- 
ular value, given the current observations and assuming 
that the reduction in flux due to Doppler smearing is cor- 
rected perfectly. To illuminate the effect of the Doppler 
smearing we calculate the average number of expected 
new discoveries akin to each of the three known DNS 
binaries. 

We have shown before (Kalogera et al. 2003) that the 
probability distribution of the expected observed number 
N^y^^ for each DNS pulsar sub-population i (B1913-I-16, 
B1534-M2, and J0737-3039) is given by: 



PiiNohs) 



Pi 



{Nobs + 1) 



(1) 



(1+A)2 (l+A)^oB.' 

where the constants Pi are a measure of how less likely it 
is to detect pulsars without acceleration searches relative 
to with acceleration searches. For each sub-population, 
the mean values of A^obs can be calculated and we find 
them to be: 

(A^'obs) 1913 = 0.9, (Arobs)l534 = 1.2, (Arobs)o737 = l-Q- 

(2) 

As expected, it is evident that the discovery of DNS pul- 
sars in tight binaries like J0737-3039 would be most fa- 
vored with acceleration searches. 

Following Kalogera et al. (2003), we can also calculate 
the combined probability distribution of the expected 
number of DNS pulsars that can be detected with PMB 
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Fig. 2. — Probability density function of the predicted number 
of observed DNS binary systems A^obs for the PMB survey, for our 
reference model (model 6 in KKL). The mean value is estimated 
to be <Nobsi=4.0. 



acceleration searches in the future. The result is shown 
for our reference model in Figure |3 The average com- 
bined number is 4.0 and the discovery of up to 2 (or 4) 
DNS systems has a probability equal to ~ 30% (~ 63%). 
We conclude that, if the acceleration search can correct 
the Doppler smearing effect perfectly, then the PMB sur- 
vey could be expected to detect an average of 3-4 DNS 
pulsars with pulse profile and orbital properties similar 
to any of the three already known systems. 

5. IN-SPIRAL EVENT RATES AND CONCLUSIONS 

Estimates of DNS in-spiral rates have suffered from 
the small number of relativistic binaries known in our 
Galaxy, mainly because of the implicit assumption in 
all methods used so far that the observed sample repre- 
sents the Galactic DNS population. Here we show that 
the recent discovery of the third relativistic binary in 



the Galactic field with binary properties (pulse-profile 
and orbital characteristics) significantly different from 
those of systems previously known reveals a new sub- 
population in the Galaxy. Consequently, it leads to a 
significant (by factors of 5-7) increase of the in-spiral 
rate estimates. 

We now consider the implications of our revised rate 
estimates for the detection of these events by LIGO and 
the other upcoming gravitational wave interferometers. 
Since these instruments can detect DNS in-spirals out 
to ~ 20Mpc for initial LIGO (~ 350 Mpc for advanced 
LIGO; Finn 2001), it is necessary to extrapolate our 
Galactic event rate out to the Local Group. Using the 
standard extrapolation of our reference model out to ex- 
tragalactic distances (Phinney 1991; KNST), we find the 
most probable event rates for our reference model are 
one per 30 yr and one per 2 day, for initial and advanced 
LIGO, respectively. At the 95% confidence interval, the 
most optimistic predictions for the reference model are 
one event per 8 yr and two events per day for initial and 
advanced LIGO, respectively. However, considering the 
full set of 27 models at 95% confidence interval indicates 
that the respective rates can reach up to one event per 3 
yr and five events per day, respectively. These results are 
quite encouraging, since, for initial LIGO in particular, 
this is the first time that DNS coalescence rate estimates 
are within an astrophysically relevant regime. Within a 
few years of LIGO operations, it should be possible to 
directly test these predictions and, in turn, place better 
constraints on the properties of binary radio pulsars and 
the cosmic population and evolution of DNS binaries. 

We also find that there is a significant probability (in 
excess of ~80%) that when acceleration searches of the 
PMB survey are completed more than 2 binary pulsars 
could be detected. The increase of the observed sample is 
very important for the reduction of the uncertainties as- 
sociated with the in-spiral rate estimates. We note, how- 
ever, that the discovery of new systems that are similar 
to the three already known does not necessarily imply 
a significant increase in the rate estimates. Significant 
changes are expected in the case that new systems with 
pulse profiles or binary properties significantly different 
are discovered, as it is such systems that will reveal a 
new DNS sub-population in the Galaxy. 
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NOTE FROM THE ERRATUM 

In our original paper, we calculated the likely size of the Galactic DNS population in two stages. First, we simulated 
the DNS distribution in the Galaxy. At this stage, in addition to storing the spatial properties and luminosities of 
the model pulsars, we also computed their expected dispersion measures and pulse scatter-broadening times using a 
model for the electron density distribution. This information was stored to an intermediate file for subsequent analysis 
by our simulation code in which the simulated population was "searched" using detailed models of the various pulsar 
surveys. Unfortunately, while the scatter-broadening times created in the first part of the calulation were saved in 
units of milliseconds, the survey simulation code assumed them to be in seconds. This error led to an under-prediction 
of the number of pulsars in the model observed samples, which in turn resulted in an overestimate of the true number, 
and hence in-spiral rate, of DNS binaries in the Galaxy. 

Wc have now repeated these calculations using the correct unit conversion and find the Galactic DNS in-spiral 
rate IZ to be 83.0lgg^'{^ Myr~^ at a 95% confidence interval (model 6), a factor of ~2.2 smaller than in the original 
paper (see Fig. 1). The corresponding detection rates for the initial and advanced LIGO detectors are: T^det.ini = 
(34.8^27:7) X 10^'^ yr^^ and 7?.dct,adv = 186.8ti48 y Y^""^) respectively. For all models we consider, the estimated DNS 
in-spiral rates range between 4 — 224 Myr~^. 

However, we note that the increase rate factor (IRF) due to the discovery of J0737-3039 remains essentially unaffected 
as the imit error cancels out, since the IRF is a measure of relative changes in rate. The IRF is found to be in the 
range 5 — 7 for all models under consideration. The revised results are summarized in Table 1. 

This error also propagates through to our predictions for future DNS discoveries in the Parkcs multibeam survey. 
Repeating the analysis in §4 of our original paper, we find that the average number of DNS with properties similar to 
those currently known to be detected using full acceleration search processing is about four. The revised probability 
distribution of expected detections is shown in Fig. 2. 
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